function  plot_gp()

load chirps.mat


s = 0.1;
th = 1;
l = 0.1;

X = chirps(1:15, 1);
t = chirps(1:15, 2);

range_start = 12;
range_end = 22;
step_size = 0.1;

range = range_start:step_size:range_end;
results = zeros(size(range, 2), 1);

for i = 1:size(range, 2)
    results(i) = gp(X, t, range(i), s, th, l).F;
end

plot(range, results)
hold on;
plot(X, t, 'g*')
print ('-dpdf','l0.1.pdf')
